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ABSTRACT 



Aims. We analyze the emission in the 0.3 - 30 GeV energy range of Gamma-Ray Bursts detected with the Fermi Gamma-ray Space 
Telescope. We concentrate on bursts that were previously only detected with the Gamma-Ray Burst Monitor in the keV energy range. 
These bursts will then be compared to the bursts that were individually detected with the Large Area Telescope at higher energies. 
Methods. To estimate the emission of faint GRBs we use non-standard analysis methods and sum over many GRBs to find an average 
signal which is significantly above background level. We use a subsample of 99 GRBs listed in the Burst Catalog from the first two 
years of observation. 

Results. Although mostly not individually detectable, the bursts not detected by the Large Area Telescope on average emit a significant 
flux in the energy range from 0.3 GeV to 30 GeV, but their cumulative energy fluence is only 8% of that of all GRBs. Likewise, the 
GeV-to-MeV flux ratio is less and the GeV-band spectra are softer. We confirm that the GeV-band emission lasts much longer than 
the emission found in the keV energy range. The average allsky energy flux from GRBs in the GeV band is 6.4 • 10~ 4 erg cirr 2 yr~' or 
only ~ 4% of the energy flux of cosmic rays above the ankle at 10 18 6 eV. 

Key words. Gamma-ray burst: general - Methods: statistical - Surveys 



1. Introduction 

Since its launch in June 2008 the Fermi Gamma-ray Space 
Telescope has broadened our understanding of Gamma-Ray 
Bursts (GRBs). The two main instruments on board of the satel- 



lite ar e the Gamma-Ray Burst Monitor (GBM, Meegan et al 



(2009)) and the Large Area Telescope (LAT, lAtwood et al 



(2009)). Together they are capable of observing GRBs over 7 
decades of energy. Especially the LAT that covers the high- 
energy region from 30 MeV to 300 GeV has the potential to pro- 
vide new insight into the underlying physics of GRBs. However 
only a small fraction of the GRBs d etected with the G BM were 
individually detected with the LAT dGehrels et al.ll2009h . 

Whereas the keV-MeV emission may well be quasi-thermal 
emission, i.e. of photospheric orig in (Eichle r & Levinsonl l2000: 
iRvde & Pe'"erll2009l: iPe'Er et al.ll2012l) . the GeV-band emis- 
sion is indicative of particle acceleration to very high energies 
and can potentially probe whether or not GRBs are powerful 
enou gh to provide a significant part of ultrah igh-energy cosmic 
rays (iKotera & Olintol201 UlEichler & Pohll201 ll) . Recently, the 
GeV-band photon output of GRBs was estimated on the basis 
of the relatively few GRBs, that are individually detected with 
Fermi-LAT, and was found to be small in comparison with the 
source power needed to s ustain cosmic rays above the ankle in 
the spectrum, at 10 18 6 eV (lEichler e t al. 2010). One source of un- 
certainty in this statement is the unknown level of high-energy 
photo n output of the many GRBs not individually detected with 
LAT (IWaxmanll2010l) . 

Here we re-analyze data of the Fermi-LAT detector with a 
view to infer the GeV-band high-energy emission from GRBs. In 
earlier studies the focus had been placed on studying each burst 
individually (e.g. | Abdo et al.ll201lt lAckermann et al.ll2012d : 
IZheng et al.l 120 1 2allbh . However, by analyzing many GRBs to- 
gether one is able to obtain more precise results for the emis- 



sion of GRBs. Therefore, we placed the focus on bursts that do 
not show a significant signal in the GeV range when analyzed 
individually. A subsample of 99 GRBs has been defined that 
were in principle detectable with Fermi-LAT and occurred in 
regions of low background emission coming from other galac- 
tic and cosmic sources. For each burst the spectrum of expected 
background photons and an effective area is estimated. Together 
with the spectrum of observed photons we thus determine the 
fluence from these bursts. 



2. Framework 

Most of the GRBs detected by the GBM do not trigger the LAT. 
Therefore it is expected that most of the GRBs do not show a sig- 
nificant signal above background. The bursts that did not trigger 
the LAT will be referred to as GBM-detected bursts. We deter- 
mine their emission by counting LAT-detected photons within a 
certain time interval and a certain solid-angle element, hence- 
forth referred to as the Area of Integration (AOI). The choice 
of the time window of observation an d the AOI is impo rtant, 
because GeV emission can be delayed ( Ab do et al .20103) and 
the GRB position is typically not well determined by the GBM. 
For long observation times and a large AOI the signal might not 
be distinguishable from the background while for short observa- 
tions and a small AOI the emission could simply be missed. A 
few general considerations are therefore in order. 



2.1. Energy Intervals 

In general the LAT is designe d to detect p hotons in an energy 
range of 30 MeV to 300 GeV (lRanddl2009h . The analysis could 
in principle be performed for the entire energy range, but the 
sensitivity of the detector varies with energy. The effective area 
of the LAT has been derived through Monte Carlo simulations 
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and verified with flight data (Ackerma nn et al.ll2012al) . and it is 
is very low for energies around 100 MeV in comparison to GeV 
energies. Thus at low energies the number of detected photons 
is low, but they represent a large portion of the flux because the 
actua l photon flux is normally a decreasing function of the en- 
ergy dAckermann et al.ll2012bl) . Additionally, the reconstruction 
of the photon arrival direction is less precise for lower energies, 
which further reduces the signal-to-background ratio. As a result 
the detection significance for energies below < 300 MeV should 
be very low and the flux determination uncertain. We therefore 
analyze the energy interval from 300 MeV to 30 GeV only. 

2.2. Time Intervals 

Recent studies of the GeV-band emission from GRBs suggest a 
significant production of highly energetic ph otons long after the 
prompt emission at keV-MeV energies (e.g. lAbdo et ai1 l2009: 
Rubtsov e t alj|2012h . mandating that the GRBs be monitored for 
long time intervals after the prompt emission phase. The dura- 
tion of GRBs is characterized by 7*90, the mid time in which 90% 
of the fluence is observed in the BATSE energy range (50 keV 
to 300 keV). The duration of the GBM emission phase and the 
LAT emission phase is likely to be correlated and therefore it 
is reasonable to choose time intervals as multiples of 790 rather 
than absolute time periods. 

2.3. Area of Integration 

From the GBM and other observatories the position of every 
burst is known with a precision of up to 1 while data from other 
observatories or a combined analysis may permit a localization 
in the arcsecond range. A comparison of such localizations in- 
dicates a systematic error in the positioning of the GRBs with 
GBM. The best fit for the systematic error is the combination 
of 2 Gaussians with dispersion 2.6° with 72% weight and 10.4° 
with 28% weight, to be added in quadrature to the statistical error 
(Paciesas et al.ll2012h . The expected distribution of GBM burst 
localization errors are shown in Fig. [1] For all other localization 
sources we neglect systematic errors. 



0.25 r 



0.20- 




Fig. 1. The expected GBM localization accuracy for 1° (dashed), 
5° (fine dashed) and 10° (continous) statistical error radius. A 
normal distribution is assumed for the statistical error. 

The imperfect angular reconstruction of the LAT, charac- 
terized by the point-spread function (PSF), may lead to a loss 
of events, even if the actual GRB was located within the AOI. 



Therefore the radius of the AOI will be the statistical error added 
in quadrature with 2° (PSF, 68% containment for 300 MeV) and 
3° (systematic error). In the case of GRB positioning with other 
detectors than GBM the systematic uncertainty is set to 0°. As 
a figure of merit roughly ~ 2/3 of all observed GeV photons of 
the GBM-detected bursts are expected to lie within the angular 
radius 



caoi 
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PSF 



(1) 



As we sum the number of gamma-ray events over many 
GRBs, the global photon-detection efficiency is approximately 
2/3. 

3. Sample Definition 

For the analysis we use a sample of GRBs detected by the 
GBM from the beginning of normal science operation of the 
Fermi satellite on A ugust 04, 2008 until July 9, 2010 as listed in 
Paciesas et alJ d2012h . In this time interval 472 bursts triggered 
the GBM. Some of these GRBs have to be excluded because they 
did not lie in the Field of View (Fo V) of the LAT or were not ob- 
served during normal science operation of Fermi. It is also useful 
to exclude certain GRBs that are expected to have a low signal 
to background ratio. After applying all cuts 99 bursts remain. 14 
of them have been previously detected by the LAT. 

3.1. LAT Field of View 

The GBM can detect GRBs in almost every direction. On the 
other hand the LAT has a FoV of ~ 2.4 sr at 1 GeV after all 
analysis cuts for background rejection have been made (Atwood 
et al. l2009h . This means that only some of the GRBs detected 
by the GBM are actually in the FoV of the LAT and therefore 
detectable by it. The sensitivity in terms of the effective area is 
a decreasing function of the angle 6 between the photon arrival 
direction and the LAT boresight. Photons arriving at 9 > 75° 
likely escape detection. It is therefore useful to concentrate only 
on GRBs where the angle between the GRB position and the 
LAT boresight during the observation was < 70°. 

3.2. Data Quality 

Some GRBs have to be rejected because of the quality of the 
data even if they are in the FoV of the LAT. As recommended by 
the LAT team, Fermi is required to be outside the South Atlantic 
Anomaly (SAA) and in normal science data-taking mode dur- 
ing the observation time. Additionally, an Earth-relative zenith- 
angle cut of 100° is imposed to prevent the spill-in of emission 
from the Earth's limb. All bursts are required to fulfill this cri- 
terion at least in a time interval of 10 x Tyo starting from the 
beginning of the T90 time. 

3.3. Background Emission 

To avoid a large number of background photons it is reasonable 
to exclude certain regions of sky that are known to have a par- 
ticularly high background intensity in the high-energy region. 
The dominant background sources are diffuse emission from the 
galactic plane as well as bright gamma-ray sources like pulsars. 
It has been found that the photon intensity at GeV energies in- 
creases by rou ghly two orders of magni tude when observing the 
galactic plane ( Ackerma nn et al.l20 1 2bl) . Regions with a galactic 
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latitude b of \b\ < 5° + o-aoi will be excluded, where q-aoi is the 
radius of the AOI specified in equation Q] It is also useful to ex- 
clude regions around very bright gamma-ray sources. The most 
luminous ones are the Blazar 3C454.3, the Vela Pulsar, the Crab 
Pulsar, the Geminga Pulsar, and the pulsars PSR J 1709-4429 and 
PSR J202 1+4026. The angular separation to these sources is re- 
quired to be at least 1° + q-aoi- Except for the Blazar 3C454.3 all 
these sources lie within the galactic plane and have already been 
excluded. 



3.4. Localization Error 

Finally it is also useful to exclude GRBs with a high localiza- 
tion error because a large AOI with correspondingly high total 
background would have to be considered. By requiring that the 
statistical error radius be smaller than 5° the number of samples 
is decreased by a factor of ~ 2/3. There is a correlation between 
the statistical error radius and the fluence. Therefore, by apply- 
ing this cut on the error angle mostly faint bursts are rejected. 



4. Method 

4. 7 . Modeling the Background 

The Fermi Science Tools are used for the analysis of the bursts. 
Although for the duration of a single GRB the background emis- 
sion is almost negligible, it becomes crucial when summing over 
many GRBs. The main sources of background radiation are the 
galactic and extragalactic diffuse emission, misclassified events, 
and individual sources like pulsars. 

Using binned likelihood analysis we construct a model of 
instrumental background and the gamma-ray sky by fitting data 
obtained in a time interval of 3 x 10 6 s until 1 hour before the 
trigger time of each burst. The pointing history of the spacecraft 
during the emission phase of the GRB then permits a prediction 
of the expected background for each burst. The event files and 
the spacecraft files are available at the Fermi mission website Q. 
The LAT has different instrument response functions (IRFs). The 
IRFs differ especially in the efficiency in terms of the effective 
area, the PSF and the energy dispersion. As recommended by the 
LAT team the P7SOURCE_V6 class is used together with the 
corresponding galactic diffuse model (gal_2yearp7v6_v0) and 
isotropic spectral template (iso_p7v6source). It offers very high 
data quality with low background from misclassified photons 
and a fairly large effective area. 

The Region of Interest (ROI) for modeling the background 
emission is a circle around the burst location with a radius of 
22°; the energy range is 300 MeV to 30 GeV. For the analysis 
only time intervals are used during which the LAT was not in 
the SAA, in normal science-operation mode, and in the normal 
range of rocking angles of not more than 52°. Additionally the 
ROI is required to not overlap with the Earth's limb. The point 
and extended sources are listed in the Fermi LAT Se cond Source 
Catal og which is based on two years of observations (Nola n et alj 
l20l2h . As recommended by the LAT team ten energy bins per 
decade and angular pixels of size 0.2 degrees are used. A first 
fit is performed with the DRMNGB optimizer and the result is 
then used for a second fit with the NEWMINUIT optimizer to 
get more precise results. 



4.2. Analysis of the Emission Phase 

For the analysis of the emission phase the P7SOURCE_V6 class 
is also used. The number of expected background photons is 
calculated by computing an exposure map for the AOI. In this 
case the exposure map describes how each source at the sky 
contributes to the observed photon counts inside the AOI. It is 
defined as: 



e(E,p) 



j ■ 

J AOI 



dt dp oh!l R(E, p, p obs ) . 



(2) 



E is the energy, p and p b s the true and observed direction 
and R is the response derived from the IRFs. This can be used 
to estimate the background photons inside the AOI from the 
background model. Energy dispersion has been neglected in this 
step. To get an estimate on the fluence of the GRB the exposure 
has to be estimated. The possible GRB positions are distributed 
throughout the AOI. The best estimate for the exposure of the 
GRB is the integral of the exposure convolved with the distribu- 
tion function <p of the expected localization errors (see for exam- 
ple Fig. [!)■ 



e(E) 



Jsky 



dp e(E, p) 



<KSp) 
2n sin(<5/)) 



(3) 



In this way the error of the burst localization and the fact that 
some GRBs actually lie outside the AOI is automatically taken 
into account. 



5. Results 

It is useful to first look at the raw photon counts in comparison 
to the number of expected photons for GRBs not detected by the 
LAT. The results are shown in Table Q] 



^start [Tgo] 


r s top[79o] 


GRBs 


"obs 


A 


-5 


-2.5 


75 


4 


2.35 


-2.5 





79 


1 


2.90 





2.5 


85 


13 


3.53 


2.5 


5 


85 


7 


3.68 


5 


7.5 


85 


8 


3.83 


7.5 


10 


85 


11 


3.93 


10 


12.5 


79 


5 


3.26 


12.5 


15 


76 


5 


2.99 


15 


17.5 


76 


6 


3.02 


17.5 


20 


73 


3 


2.92 


20 


22.5 


71 


8 


2.36 


22.5 


25 


67 


4 


2.17 


25 


27.5 


64 


3 


1.92 


27.5 


30 


64 


5 


1.88 


30 


32.5 


63 


1 


1.78 


32.5 


35 


61 


1 


1.57 





10 


85 


39 


14.97 



Table 1. The number of observed photons (« bs) and the ex- 
pected background {A) for the GBM-detected bursts for the dif- 
ferent time intervals. The times are given in multiples of the 790 
time of each burst and relative to the start of the Tw time inter- 
val. For times outside the 0-10 Tgg time interval some GRBs 
have to be excluded on account of cuts that were originally only 
applied in this time interval. 



http : //heasarc . gsf c . nasa . gov/ 
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As we can see there is a significant emission in the time in- 
terval from to 10 times the Tyo time. In this time A = 14.97 
photons are expected and n — 39 observed. The chance proba- 
bility for a background fluctuation is 1 .7x 10~ 7 , based on Poisson 
statistics 



ZA n 
~l e ~ A 



(4) 



This is a significant signal when summed over all GRBs but 
clearly not significant for a single burst for which the number of 
observed counts is typically smaller than 1 and the number of 
observed photons mostly or 1 . For the 1 1 observed short bursts 
(T90 < 2 s) there was 1 observed photon and 0.06 expected back- 
ground photons which is too low a signal to claim a detection. 
For the remaining 74 long bursts (790 > 2 s) there were 38 ob- 
served photons and 14.91 expected from background. Note that 
before the GRB trigger we observe 5 events and expect 5.25. 
Later than 25 Tgo after trigger we observe 10 events and expect 
7.15, indicating that there is little, if any, emission at that late 
stage, and no evidence for any activity before the GBM trigger. 

5.1. Detection Significance 

One can search for significant signals of single GRBs in these 
data. For this purpose the number of observed and predicted pho- 
tons are compared for the 4 time intervals between to 10 Tgo as 
well as the total time interval. As a cross-check the LAT-detected 
bursts will also be tested with the same method. 99 bursts are 
tested for 5 different time intervals and therefore we have at 
most 495 trials, because the time intervals are not all indepen- 
dent. Therefore the probability that the background emission, A, 
reaches or exceeds the observed emission, « D b s , by normal fluc- 
tuations should be smaller than 10~ 4 for a 5% post-trials chance 
probability. All 14 bursts previously detected by the LAT show 
a significant signal in at least one time interval. This is not sur- 
prising as these bursts actually triggered the LAT and also have 
high-precision localizations either by the LAT or other observa- 
tories. However there were 2 other bursts with significant emis- 
sion: GRB 100207B and GRB 08 1009 A. The latter one has also 
been reported as a candidate for a LAT-detected burst with the 
same method but with a wider energy range (E > 100 MeV) an d 
longer observations times (f D b s = 1500 s) (Rubtsov e t alJl2012h . 
However it has not been added to the LAT Burst catalog so far. 
The results of this analysis are shown in Table [2] 



GRB 


Time [r 90 ] 


"obs 


A 


P(n > n obfi , A) 


08 1009 A 


5-7.5 


3 


0.0366 


7.93 x 10~ 6 


100207B 


0-2.5 


2 


0.0123 


7.54 x 10- 5 



Table 2. We show results for the two GRBs not listed in the 
LAT Burst catalog that show significant emission in the analysis. 
As before, n j, s is the number of observed photons and A the 
expected number of photons from background. 



5.2. Photon Fluence 

The average photon fluence of GRBs can be calculated by sep- 
arately summing the exposure and the observed photon and 
expected background counts in every energy bin. To estimate 
the cumulative exposure over the entire energy range between 



300 MeV and 30 GeV we assume the emission follows a power- 
law with index of -2.3. As we shall see below, this index is ap- 
propriate for LAT-detected GRBs, whereas the average emission 
of GBM-detected bursts is better described with a power-law in- 
dex of -3, in which case the photon fluence would be underesti- 
mated by about 10%. The 68%-confidence statistical uncertainty 
of the fluence has been calculated using an incomplete gamma 
function, i.e. the Bayesian inversion of a Poisson distribution for 
uniform prior, and is considerably larger than the systematic un- 
certainty in estimating the fluence. The analysis has been done 
for the GBM- and LAT-detected bursts separately. The result is 
shown in Tableland graphically displayed in Figure|2] 



Time [T 90 ] 


Fluence [y 


cm~ 2 ] 




GBM 


LAT 


(-5) - (-2.5) 


°- 93 -!:3o x 10 ~ 5 


-0.64+ 347 x 10~ 5 


(-2.5) -0 


-10.29+1^ x 10~ 6 

— 3.4U 


5.34+^™ x 10- 5 


0-2.5 


4.86+|fj x 10- 5 


66.03+ 4 ;} 2 x 10- 4 


2.5-5 


1.67+J-3' x 10- 5 


9.63+Jj 2 x 10- 4 


5-7.5 


2.07+jj® x 10~ 5 


5.75+\f g x 10- 4 


7.5-10 


3.47!}$ x 10~ 5 


15.34+I;* 2 x 10- 5 


10-12.5 


0.92+}$ x 10- 5 


33.64+J;® x lO" 5 


12.5-15 


1.08+j ^ x 10-5 


18.00+| 44 x 10- 5 


15-17.5 


i.6o+J:S x 10-5 


13.71+1;™ x 10~ 5 


17.5-20 


0.05:3$ x 10-5 


11.75!»Jxl0- 5 


20 - 22.5 


3.18+}$ x 10- 5 


14.42+J;^ x 10~ 5 


22.5-25 


1.09 + | 11 x 10- 5 


16.49^ x 10- 5 


25-27.5 


0.67+1* x 10~ 5 


1.96+|g x 10- 5 


27.5-30 


L95 -L55 x 10_5 


1.95+||° x 10~ 5 


30-32.5 


-5.00+j$ x 10- 6 


4.48+^ x 10- 5 


32.5-35 


-3.77+^x10"' 


2.06+|;« x lO" 5 



Table 3. We show the photon fluences for different time inter- 
vals around the start of the Tgo time interval. The GBM fluence 
is the weighted average fluence of the 85 bursts not detected by 
the LAT. The LAT fluence is the weighted average fluence of 
the 14 bursts detected by the LAT, which are individually listed 
in Table [4] The LAT fluence is dominated by the three excep- 
tionally bright bursts GRB 0809 16C, GRB 090902B and GRB 
090926A. 



Both GBM- and LAT-detected bursts show GeV-band emis- 
sion long after the bulk of emission in the BATSE energy range 
(50 keV to 300 keV), which arises during the Tgo time. Whereas 
the high-energy photon fluence in the Tgo time interval from the 
GBM-detected bursts is 2.38+} 2 j X 1CF 5 y crrT 2 , for the time 
interval from to 25 times T^o we find an average photon flu- 
ence from the GBM-detected GRBs of 19.98+^ x 10~ 5 y cm 2 . 
Given that there could still be emission at later times which is 
just too low to be detected, the observed fluence during the Tgo 



time is at most 1 1 .91+5'5Z% of the total GeV-band emission. One 
has to keep in mind that this is the average fluence weighted 
with the effective area of each burst. Since the Fermi satellite 
occasionally reorients for a very bright burst, our method may 
favor bright bursts especially for later times. 



4 



J. Lange & M. Pohl: GeV-band Emission from GRBs 







I 

1 44 







10 20 30 

Time [r 90 ] 

Fig. 2. Average fluences for different time intervals in units of 
Tgo. Shown in black (grey) are the average fluences for the 
GBM-detected (LAT-detected) bursts. The fluences for the LAT 
detected bursts were scaled by a factor of 10~ 2 to fit on the same 
plot. 



The photon fluence has also been separately calculated and 
listed in Table |4] for all LAT-detected bursts using the data from 
to 10 times the Tgo time interval. The results thus derived are 
comparable to thos e previously found with standard methods. 
The Burst Catalog (Paciesasetal .2012) can be used to compare 
the fluence in the high-energy region from 300 MeV to 30 GeV 
to that in the BATSE range from 50 keV to 300 keV obtained 
with the GBM. To be noted from Table [4] is that the bursts not 
detected by the LAT are fainter in the high-energy region than 
the LAT-detected GRBs when normalized to the same keV-band 
fluence. 



GRB 


^batselT crrT 2 ] 


F [300 MeV,30 GeV] [j Cm 2 ] 


080825C 


114.2 ±0.2 


1.81 + };f° x 10~ 3 


0809 16C 


6.93 ± 0.2 


x 10~ 3 


081024B 


0.16 ±0.03 


12.93^° x 10- 4 


090217 


4.8 ±0.1 


2.12!^ X10- 3 


090328 


7.8 ±0.2 


3.59 + |; 2 ? x 10~ 3 


090510 


1.11 ±0.05 


11.99+}$! x 10~ 3 


090626 


15.0 ±0.3 


2.19 + |; 3 «x 10- 3 


090902B 


29.5 ± 0.3 


35.87+ 3 ;f» x 10~ 3 


090926A 


40.1 ± 0.4 


26.64 +3 24 x 10~ 3 


091003 


17.1 ± 0.3 


9 "-I:o5 x 10-4 


091031 


3.1 ±0.2 


2.55!};« x 10- 3 


100225 A 


1.9 ±0.1 


0.95 + _' f 6 x 10- 3 


100325 A 


2.6 ±0.1 


4.39 + ^ x 10"* 


100414A 


10.1 ± 0.4 


6.83 +3 ;° 3 x 10~ 3 


All GBM Bursts 


3.38 ±0.01 


19.98 +4 ^ 2 x 10- 5 



Table 4. Comparison of the fluences from several LAT-detected 
bursts and the average fluence for the 85 GBM-detected bursts 
with the fluences in the BA TSE range taken from the GRB 
Catalog dPaciesas et alJl2012D and derived from a fit to a Band 
function. 



5.3. Average Flux from GRBs 

We now determine the average emission from GRBs in the en- 
ergy range from 300 MeV to 30 GeV. GRBs not detected by 
the LAT contribute with a fluence of 17+ 3 - 7 x 10~ 3 y cm' 2 , 
whereas LAT-detected GRBs account for a fluence of 120+Zg x 
10~ 3 y cirT 2 . Hence the total observed fluence in the range from 
300 MeV to 30 GeV is 137 + j^ x 10~ 3 y cirT 2 . Since 50% of the 
observed fluence comes from only 2 GRBs (GRB 090902B and 
GRB 090926A) the statistical error on this estimate should be 
roughly 0.5/ V2 * 35%. 

To estimate the average allsky flux of GRB-produced GeV- 
band gamma rays one has to consider all cuts applied in the anal- 
ysis. First, some bursts were dropped on account of data quality 
as discussed in section [3~2l This affected 27% of all available 
bursts giving a weight factor of 1.37. Then, the 5° cut on the 
statistical error radius of the burst location is compensated by 
an additional factor. The bursts excluded in this step account for 
3% of the overall fluence observed in the GBM energy range. 
This gives a factor of 1.03, assuming that the fluences in the 
GBM energy range and in the LAT energy range are propor- 
tional. (We find less GeV emission than that, and therefore the 
true correction factor is somewhere between 1 and 1 .03). Finally, 
the limitation to bursts that occurred in the LAT FoV (0 < 70°) 
and were not too close to the galactic plane and the Blazar 
3C454.3 is compensated with a factor of 3.46. We combine the 
three correction factors to derive the total efficiency factor as 
4.88 = 3.46 • 1.03 • 1.37. The total high-energy fluence over the 
entire sky should therefore be (6.69 + 2.37) x 10 _I y crrT 2 . The 
effective observation time from August 04, 2008 until July 9, 
2010 was 5.083 x 10 7 s, and hence the total averaged flux from 
all bursts over the entire sky in the energy range from 300 MeV 
to 30 GeV is (13.2 ± 4.7) x 10~ 9 y cirT 2 s _1 . 

We have also determined the energy fluence from GRBs. 
GRBs not detected by the LAT contribute with a fluence of 
1.67 x 10" 5 ergcm~ 2 , whereas LAT-detected GRBs account for 
a fluence of 1.93 x 10~ 4 erg cirT 2 . The average photon energies 
are ~ 1 GeV for LAT-detected GRBs and ~ 0.6 GeV for GBM- 
detected bursts, estimated as the ratio of the energy fluence and 
the photon fluence. If the emission spectra were power laws, then 
the spectral indices, s, would be slat — 2.3 and sgbm — 3. 

The average allsky energy flux is about 6.4 • 
10~ 4 ergcm~ 2 yr^' , slightly larger than the estimate of 
lEichler et al. (2010), but still compatible with it considering the 
uncertainties. 

6. Summary and Discussion 

In this study of high-energy emission from GRBs in the range 
from 300 MeV to 30 GeV the focus has been placed on GRBs 
not individually detected by the LAT and referred to as GBM de- 
tected. We defined a sample of 85 GRBs listed in the GBM cata- 
log, 74 of which can be regarded as long bursts (Tqo > 2 s) while 
the remaining 11 GRBs are short (790 < 2 s) bursts. We find 
significant emission above background for the complete sam- 
ple. For the long bursts the emission is clearly visible while 
for the short bursts alone the results are statistically inconclu- 
sive. Moreover, the GeV-band emission lasts considerably longer 
than the 790 time of the keV-MeV energy range, in fact at least 
to 10 times T^. A similar conclusion was previously reporte d 
for individual GRBs like GRB 80916C (lAbdo et alJ 12 009). 
GRB 081024B ( lAbdoetalJ2010bl) or GRB 940217 (lHurlevetal.l 
119941) . Altogether only (12 + 7)% of the total photon fluence in 
the range from 300 MeV to 30 GeV is emitted during Too. Since 
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this extended emission was observed when summing over many 
GRBs, we cannot discriminate between continous emission and 
a sequence of flares. Given that the number of GRBs in the sam- 
ple is larger than the number of observed photons, it is also un- 
clear whether this delayed emission is a general feature or the 
product of a subset of GRBs. 

The ratio of the fluence in the high-energy region to the flu- 
ence in the BATSE energy region is lower (6 • 10~ 5 ) for the bursts 
not detected by the LAT than for LAT-detected bursts (4.7 • 10" 4 ). 
The estimated spectra in the GeV band are softer for GBM- 
detected bursts (with equivalent photon index s ^ 3) than for 
LAT-detected bursts (equivalent photon index s 2.3). There 
is no indication that GRBs produce particle spectra with typi- 
cal indices s ^ 2 which are often assumed in source models of 
ultra-high-energy cosmic rays. Likewise, there is no evidence for 
a population of GRBs that efficiently accelerate particles to high 
energies but individually emit too weakly in the GeV band for a 
detection with Fenni-LAI. 

Altogether the bursts not detected by the LAT contribute 
roughly 14% of all GRB-produced photons, and 8% of the emit- 
ted energy, in the energy range from 300 MeV to 30 GeV. 
Finally, we find that the average allsky gamma-ray flux com- 
ing from GRBs in this energy range is (13.16 ± 4.65) x 
10~ 9 y cirT 2 sT l . The average allsky energy flux from GRBs in 
the GeV band is only ~ 4% of the energy flux of cosmic rays 
above the ankle at 10 18 6 eV. 
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